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摘 要 : 全 球 城市 化 进程 加 快 引起 的 干旱 区 绿洲 城市 热岛 效应 变化 及 其 生态 环境 问题 已 成 为 当前 
国内 外 城市 气候 、 环 境 和 生态 等 领域 的 研究 热点 之 一 。 基 于 遥感 热 红 外 通道 反 演 的 地 表 温 度 
(Land surface temperature ,LST) 是 开展 城市 地 表 热 岛 (Surface urban heat island,SUHI) 效 应 监测 研 
究 的 关键 参数 。 然 而 ,受热 红外 遥感 数据 "时空 矛盾 "的 制约 ,目前 在 单一 星 载 卫星 传 感 器 下 尚 不 
存在 同时 具有 高 时 间 分 辨 率 和 高 空间 分 辨 率 的 热 红 外 遥感 数据 源 , 因 而 制约 了 在 干旱 区 绿洲 城市 


范围 内 开展 高 精度 地 表 热 岛 效应 监测 研究 。 针 对 上 述 问题 ,以 干旱 区 典型 绿洲 城市 


乌鲁木齐 


市 为 研究 区 域 , 以 Landsat 系列 影像 和 MODIS 地 表 温 度 产 品 为 基础 数据 源 , 基 于 FSDAF (Flexible 
spatiotemporal data fusion method ) 时空 融 合 模型 分 析 了 2001—2018 年 乌鲁木齐 市 在 城市 扩张 背景 
下 夏季 地 表 热 岛 效 应 的 时 空 变化 特征 以 及 夏季 LST 与 城市 地 表 参 数 之 间 的 关系 。 研 究 结果 表明 : 
(D 乌鲁木齐 市 夏季 热岛 强度 (SUHI intensity,SUHII) 在 不 同 的 郊区 范围 内 均 呈 现 增加 的 趋势 。 在 
较 小 的 郊区 范围 内 ,SUHIL 从 2001 年 的 1.24 % 增 加 到 了 2018 年 的 2.83 %C; 在 较 大 的 郊区 范围 内 ， 
SUHID 从 2001 年 的 1.44% 增 加 到 了 2018 年 的 2.88%C;(2) 在 研究 区 各 土地 利用 类 型 中 , 裸 地 的 夏 


季 LST 最 高 ,水 体 最 低 ;(3) 研究 区 地 表 反 照 率 和 不 透水 
而 植被 指数 与 植被 覆盖 度 则 与 LST 呈 负 相 关 关 系 ;(4) 在 干旱 区 绿洲 城市 ,城区 内 部 植被 


的 增加 与 城市 夏季 LST 升 高 呈正 相关 ， 
积 的 增 


加 有 助 于 缓解 城市 热岛 效应 ,而 仅 郊区 植被 的 增加 则 会 导致 SUHI 的 进一步 增强 。 
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近年 来 ,热岛 效应 随 着 全 球 城市 化 进程 的 不 断 
加 快 ,目前 已 成 为 一 个 影响 城市 居民 身心 健康 , 城 
市 空气 质量 ,水 资源 消耗 和 城市 生态 系统 可 持续 发 
展 的 重要 生态 环境 问题 "。 相 对 于 利用 传统 气象 站 
观测 数据 和 数值 模拟 开展 的 城市 大 气 热岛 监测 , 星 
载 和 机 载 等 怕 感 热 红外 传感器 因 其 具有 罗 盖 面 广 、 
获取 数据 成 本 低 等 优势 目前 已 广泛 用 于 城市 地 表 
热岛 (Surface urban heat island, SUHI) 监测 人 研究。 
地 表 温 度 (Land surface temperature ,LST) 是 基于 热 
红外 传 感 右 反 演 得 到 的 ,可 以 实现 城市 地 表 热 驴 高 
精度 监测 的 关键 参数 。 然 而 由 于 星 载 热 红外 传 感 
顺 设 计 和 研发 技术 的 瓶颈 限制 ,目前 尚 不 存在 同时 
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满足 高 时 间 和 高 空间 分 辩 率 的 热 红 外 遥感 影像 ”。 
例如 :目前 常见 的 高 时 间 分 辨 率 MODIS 地 表 温 度 产 
品 的 空间 分 辨 率 较 低 (1 km) ,而 高 空间 分 辨 率 的 
Landsat 系列 卫星 数据 的 重 返 周期 则 长 达 16 d, 这 种 
遥感 数据 的 “时 空 矛 盾 "不仅 导致 了 高 时 空 分 辨 率 
LST 数 据 的 缺失 ,还 严重 限制 了 遥感 数据 在 区 域 地 
表 热 岛 效 应 监测 中 的 应 用 。 

目前 ,地 表 温 度 时 空 融 合 模 型 已 成 为 解决 上 述 
问题 的 有 效 手段 ,其 主要 原理 在 于 融合 高 空间 分 辩 
率 地 时 间 分 辨 率 LST 数 据 和 低空 间 分 状 率 高 时 间 分 
辨 率 LST 数 据 的 特征 ,从 而 生成 高 时 空 分 辩 率 的 
LST 数 据 。Liu 和 Weng" 采 用 时 空 自 适 应 反射 率 融 
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E RÆ. 基于 FSDAF 模型 的 干旱 区 


合 模型 (Spatio temporal adaptive reflectivity fusion 
model, STARFM) ,将 MODIS LST 24 fil ASTER LST 
数据 相 结合 ,预测 生成 的 高 时 空 分 辨 率 LST 影 像 与 
真实 影像 的 差 值 小 于 1 Co Wu 4° Landsat LST 
SCE MODIS LST 数 据 和 GEOS/SEVIRI 数 据 相 结合 
构建 了 STITFM (Spatio-temporal integrated tempera- 
ture fusion model) 模 型 ,预测 生成 了 高 时 空 分 辩 率 的 
LST 数据 。 郑 明亮 等 "基于 TsHARP (Thermal 
sHARPening) 模 型 和 STITFM 模型 预测 生成 了 30 m 
空间 分 辨 率 的 LST 数 据 , 其 均 方 根 误差 小 于 1.33 
K。 魏 然 和 单 杰 " 基 于 ESTARFM (Enhanced spatial 
and temporal adaptive reference fusion model) 模 型 和 
非 线 性 辐射 温度 分 解法 融合 生成 了 60 m 空间 分 辨 
率 的 逐日 LST 影 像 。Zhu 等 提出 了 可 以 有 效 减 小 
误差 ,从 而 预测 时 序 高 分 辨 率 LST 数 据 的 FSDAF 
(Flexible spatiotemporal data fusion method ) 模 型 。 研 
究 结 果 表 明 ,FSDAF 模 型 不 仅 运 行 效率 较 高 ,数据 
源 输 入 需求 较 小 ,预测 生成 的 LST 数 据 也 具有 较 高 
精度 。 
基于 此 ,本 研究 以 近 20 a 来 受 经 济 快速 发 展 影 
响 ,城市 化 进程 最 为 明显 的 干旱 区 典型 绿洲 城市 
一 一 乌鲁木齐 市 为 研究 区 ,基于 FSDAF 模型 将 反 演 
获取 的 Landsat LST 数 据 和 MODIS LST 数 据 相 结合 ， 
预测 生成 高 时 空 分 辨 率 的 LST 数 据 。 在 此 基础 上 ， 
开展 了 2001 一 2018 年 乌鲁木齐 市 夏季 地 表 热 岛 时 空 
变化 特征 进行 研究 ,重点 探讨 了 人 研究 区 夏季 LST 与 土 
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典型 绿洲 城市 夏季 地 表 热岛 效应 时 空 演 变 研究 


地 利用 类 型 和 城市 地 表 参 数 之 间 的 关系 ,以 期 为 今 
后 干旱 区 绿洲 城市 的 地 表 热 岛 效 应 研究 提供 参考 。 


1 研究 区 概况 


乌鲁木齐 市 位 于 我 国 西北 部 ,是 新 疆 维 吾 尔 自 
yá EX AY ET JF (86° 37’ 33" ~88°58' 2A" E, 42°45'32"N~ 
45*00'00"N) ,城市 北部 为 河流 冲击 平原 ,南部 与 天 
山中 段 北欧 相 接 ,东北 部 与 西南 部 也 均 为 山地 , JÉ 
成 三 面 环 山 的 格局 (图 1)”。 乌 鲁 木 齐 市 平均 海拔 
TE 800 m 左 右 ,属于 中 温带 大 陆 性 干旱 气候 ,年 均 降 
水 量 为 145 mm, 四季 差异 明显 ,夏季 为 每 年 6 月 中 
旬 一 9 月 中 旬 ,夏季 持续 时 间 长 ` 昼 夜 温 差 大 且 炎 热 
干燥 ,是 典型 的 干旱 区 绿洲 城市 。 

乌鲁木齐 市 作为 中 亚 最 大 的 干旱 区 绿洲 城市 
之 一 , 近 20 a 来 经 济 快速 发 展 ,GDP H 2001 年 的 
315x10* 元 增加 到 了 2018 年 的 3099x10: 元 。 同 期 常 
住人 口 数量 也 由 169.03x10!: 人 快速 增长 到 350.58x 
10: 人 。 人 口 数量 激增 原 有 自然 表面 大 量 被 人 工 不 
透水 面 所 蔡 代 ,城市 整体 热 排放 水 平 的 不 断 增加 均 
造成 了 城市 热 环境 问题 的 逐步 加 剧 。 


2 数据 来 源 和 研究 方法 


2.1 Landsat LST 数 据 获取 
本 研究 采用 的 Landsat 系列 数据 (path/row: 143/ 
29; path/row: 143/30) 的 时 间 分 辩 率 均 为 16 do H 


45°N 


44°N 


图 1 研究 区 概况 
Fig. 1 Location ofthe study area 
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表 1 本 研究 在 测试 实验 阶段 选用 的 Landsat TM、Landsat 8 和 MODI11A1 数据 
Tab. 1 Selected Landsat TM, Landsat 8 and MOD11A1 data for the testing experiment in the study 


已 知 时 刻 (5) 基 础 预测 时 刻 ( 心 ) 基 础 


(年 -月 -日 ) 时 间 (UTC) MALTE 输入 LST 影 像 

Landsat TM 2001-05-23 04:36 v 

MOD11A1 2001-05-23 04:30 v 

MOD11A1 2001-07-10 04:30 V 

Landsat TM 2001-07-10 04:36 v 
MOD11A1 2001-09-12 04:30 v 

Landsat TM 2001-09-12 04:36 v 
Landsat TM 2011-06-04 04:45 v 

MOD11A1 2011-06-04 04:20 v 

MOD11A1 2011-08-07 04:20 v 

Landsat TM 2011-08-07 04:44 v 
Landsat 8 2018-05-22 04:55 v 

MODIIAI 2018-05-22 04:25 v 

MOD11A1 2018-06-23 04:55 v 

Landsat 8 2018-06-23 04:55 v 
MOD11A1 2018-08-26 04:55 v 

Landsat 8 2018-08-26 04:55 V 
MOD11A1 2018-09-11 04:55 v 

Landsat 8 2018-09-11 04:55 v 


中 ,Landsat TM 和 Landsat 8 热 红 外 波段 的 空间 分 辩 
率 分 别 为 120 m 和 100 m。 数 据 来 源 于 U.S. Geologi- 
cal Survey 网 站 (http://earthexplorer.usgs.gov/) , 热 红 
外 波段 的 空间 分 辨 紊 已 经 在 下 载 后 重 采 样 到 了 30 
m, 研究 所 采用 的 的 Landsat 数据 清单 如 表 1 所 示 。 
在 预 处 阶段 ,数据 本 身 已 经 进行 过 了 几何 校正 ,并 
在 此 基础 上 完成 了 辐射 校正 和 大 气 校 正 。 最 后 分 
别 基 于 Landsat TM 热 红 外 第 6 波段 和 Landsat 8 热 红 


外 第 10 波段, 采用 单 窗 算法 反 演 获取 了 Landsat LST 
数据 。 有 具体 的 单 窗 算 法 反 演 的 详细 步骤 见 Qin SE”! 
的 研究 成 果 。 
2.2 MODIS 地 表 温 度数 据 预 处 理 

本 人 研究 同时 采用 当地 时 间 10:30 过 境 的 每 日 合 
成 MODIS LST 产品 MOD11A1 和 8 d 合成 MODIS 
LST 产 品 MOD11A2( 表 2) ,二 者 的 空间 分 辩 率 均 为 
1 km ,数据 来 源 于 美国 NASA 的 对 地 观测 数据 和 信 


表 2 本 研究 在 融合 实验 阶段 选用 的 MOD11A2 数 据 
Tab.2 Selected MOD11A2 data for the fusion experiment in the study 


Totg Za] C ) 8) LST HAR 


实验 年 中 天 数 /d 数据 获取 日 期 = 
2001 年 2011 年 2018 年 
融合 实验 169~176 6 H 18—25 H v v 
177~184 6 月 26 日 一 7 月 3 日 V v 
185-192 7H4—11H v M M 
193~200 7H 12—19 H v v v 
201-208 7H20—27H v V v 
209~216 7 月 28 日 一 8 月 4 日 v v v 
217-224 8H5—12H v M M 
225-232 8H 13—20 H v v v 
233-240 8H21—28 H v M M 
241-248 8 月 29 日 一 9 月 5 日 v V M 
249-256 9 月 6 日 一 9 月 13 日 v M v 
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息 系 统 (https://earthdata.nasa.gov) o MODI1A1 和 
MOD11A2 是 目前 应 用 于 大 区 域 尺 度 下 热 面 地 表 温 
度 监 测 适用 性 最 好 的 产品 ,目前 已 广泛 用 于 城市 地 
表 热 岛 监测 研究 。 首 先 , 本 研究 利用 MODIS 数据 预 
处 理工 具 MRT 将 MODIS LST 数 据 的 坐标 系统 转换 
成 与 Landsat 数据 一 致 的 UTM 坐标 系 。 其 次 ,以 
Landsat LST 数 据 为 基准 ,对 MODIS LST 数 据 进行 几 
何 校正 ,消除 不 同 传感器 数据 源 之 间 的 几何 误差 。 
2.3 FSDAF 模型 原理 

Zhu 等 中 提出 的 FSDAE 模 型 ,由 于 仅 需要 在 已 
知 时 刻 的 一 对 高 .低空 间 分辩 率 遥 感 数据 和 预测 时 
刻 的 一 景 低空 间 分 辩 率 遥感 数据 ,就 可 以 有 效 融 合 
生成 预测 时 刻 的 高 空间 变 率 影像 ,因而 目前 被 国内 外 
学 者 应 用 于 LST 数 据 的 时 空 融 合 研究 。Zhang 等" 利 
用 FSDAF 模型 预测 后 成 了 时 序 高 空间 分 辨 率 LST 
数据 集 。 杨 敏 等 王 利 用 覆盖 黑河 中 游 盘 科 灌 区 的 
ASTER LST fil MODIS LST 数 据 , 基 于 FSDAF 模 型 预 
测 生 成 了 与 ASTER 相同 空间 分 辩 率 的 LST 数据 
集 。 因 此 ,本 研究 基于 FSDAF 模型 将 高 空间 分 辩 率 
AY Landsat TM 和 Landsat 8 LST 数 据 与 低空 间 分 辨 率 
的 MOD11Al 和 MOD11A2 数 据 相 结 合 , 以 期 预测 生 
成 高 时 空 分 状 率 的 LST 数 据 。 本 研究 将 基于 FSDAF 
模型 融合 生成 的 高 时 空 分 辩 率 LST 数 据 称 为 "LST 
预测 数据 ” ,将 基于 MODIS il Landsat 卫星 反 演 生成 
的 LST 数 据 称 为 "LST 观 测 数据 ”。 

本 研究 有 2 个 重要 步 又 :测试 实验 和 融合 处 
理 。(1) 测试 实验 :为 了 验证 FSDAF 模 型 融合 LST 数 
据 的 精度 能 够 满足 研究 需要 ,本 研究 首先 利用 
2001. 2011 年 和 2018 年 夏季 开始 前 1 个 月 内 的 
Landsat LST 和 MOD11A1 数据 为 FSDAF 模 型 基础 时 
刻 (4) 的 输入 影像 ( 表 1) ,在 此 基础 上 利用 夏季 预测 
时 刻 (,) 的 MODI11A1 数据 融合 生成 夏季 不 同 预测 
时 刻 的 高 空间 分 辨 率 “LST 预 测 数据 ”。 最 后 ,利用 
基于 Landsat 影 像 反 演 得 到 的 预测 时 刻 “LST 观 测 数 
据 ”" 对 “LST 预 测 数 据 ” 的 精度 进行 验证 。FSDAF 模 
型 的 具体 处 理 步 又 如 下 所 示 : 
Lit, (x, y,) -LST,(x,, »,) + YW x ALSTœ, y) (1) 


ALST(x, Yy) 6g (x Y) + ALST() (2) 
式 中 : Lêt, (x,, y) A MT ZA, OP 28 T8] PER 
“LST 预测 数据 ; LST, Gc, V) 表示 基础 时 刻 4 的 高 
25 HAPEE LST fiii AGIS WRR AS k RTC 


的 权重 ; ALST(x,, y,) 表示 刻 到 六 时 刻 相似 像 元 LST 
的 变化 值 ; ALST(x,, y;) 表示 目标 像 元 Qs. y,) TE ti 
到 六 时 刻 的 变化 值 ; eu Goss Yi) 表示 第 1 个 高 空间 分 
辩 率 像 元 与 第 ;个 低空 间 分 辩 率 像 元 之 间 的 残 差 ; 
ALST(c) 表示 高 空间 分 辨 率 LST 数 据 在 第 类 别 上 
从 石 到 记 时 刻 的 变化 值 。 更 多 FSDAF 模型 权重 的 计 
算 步 又 和 细节 请 参考 Zhu 等 的 研究 成 果 。 
可 以 用 于 计算 残 差 的 TPS (The thin plate 

spline ) 函数 求解 方法 如 下 所 示 : 

Erigh y» y) 5m x L(x, y)* W(x;, »j) (3) 

L(x; y) - AL(x;, y)- 


二 > LST py)— DLST, (x, Vi) 
ml j=l 


(4) 


LW(x;, yj) = EX; yj) x 
HI(x,, ys) 十 五 (Xp Vi) x[I-HKx,, A 
E Xy y) - LST Qt, YLST (x, yy) (6) 


LST (x,, Paice) (7) 
式 中 : Lo, y) 表示 高 空间 分 辩 率 "LST 观 测 数据 ?和 
“LST 预测 数据 ?之 间 的 残 差 值 ; 普 表示 低空 间 分 辩 
率 影像 的 亚 像 元 数量 ; AL(x,, V) 表示 LST 在 低空 间 
分 辨 率 影像 从 4 到 志 时 刻 的 变化 值 ; LW(x,, y,) 表示 
残 差 分 配 的 权重 ; W(x. y;) 表示 LWO, y;) 的 归 一 
化 权重 ; LST)? 为 基于 时 间 变 化 预测 生成 的 时 刻 
LST 值 ; E, (x, y;) 为 时 间 的 预测 误差 ; HI y) 为 
同 质 系数 ; E, Qs, y) 表示 在 低空 间 分 辩 率 像 元 中 
的 等 误差 ; LST2 (zy y) 为 基于 优化 参数 后 的 TPS 
函数 融合 生成 的 总 时 刻 高 空间 分 辩 率 LST 影 像 ; 
Pros Qc, Vy) 为 TPS 函数 。 

本 人 研究 采用 预测 时 刻 反 演 获取 的 高 空间 分 辩 
RLST 观测 数据 ”作为 参考 影像 ,选用 决定 系数 
(R) \ 均 方 根 误差 (RMSE) 和 绝对 平均 差 (AAD ) 等 3 
个 定量 评价 指标 验证 “LST 预 测 数 据 ” 的 精度 。 如 果 
本 研究 在 测试 实验 阶段 预测 生成 的 LST 数 据 的 精度 
能 够 满足 研究 需要 ,人 研究 再 进行 融合 处 理 ,(2) 融合 
处 理 : 研 究 同样 利用 2001.2011 年 和 2018 年 夏季 开 
始 前 1 个 月 内 的 Landsat LST 4l MODI1A1 遥感 数据 
为 FSDAF 模 型 基础 时 刻 (4) 的 输入 影像 ,将 夏季 预 
测 时 刻 (4) 8 d 合 成 的 MOD11A2 数 据 替 换 原 有 预测 
时 刻 (6)MOD11A1 数据, 从 而 预测 生成 时 间 分 辩 率 
为 8d 窗 羡 人 研究 区 的 夏季 高 空间 分 辨 率 LST 数 据 集 ， 


(5) 
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并 分 别 将 该 数据 集 内 2001、2011 年 和 2018 年 的 夏 
季 LST 数 据 之 和 求 平均 ,从 而 最 终 获 取 人 研究 区 各 年 
份 夏季 平均 LST 数 据 。 
2.4 热岛 强度 (SUHI Intensity,SUHII) 计算 

研究 采用 SUHIL 和 SUHIL 等 2 个 热岛 强度 指标 
计算 研究 区 热岛 强度 。SUHIL 定 义 为 市 区 与 郊区 
(面积 占 市 区 面积 的 100% ) 夏 季 平 均 LST 温 度 差 ; 
SUHIL 定义 为 市 区 与 郊区 (面积 占 市 区 面积 的 
150% ) 夏 季 平 均 LST 温 度 差 ,以 此 验证 研究 区 不 同 
郊区 范围 对 地 表 热 岛 效应 变化 的 有 影响。 郊区 范围 
的 界定 参考 Peng 等 ”的 研究 成 果 。SUHII 的 计算 公 
式 如 下 所 示 : 

SUHII, =LST „san -LST porat 100% (8) 
SUHI, = LST svoan LST rrai 150% (9) 

2.5 城市 地 表 参 数 计 算 

为 了 进一步 探索 研究 区 夏季 IST 与 城市 地 表 参 
数 变 化 之 间 的 关系 ,研究 采用 了 增强 植被 指数 (En- 
hanced vegetation index, EVI) 植被 覆盖 度 (Fraction- 
al vegetation cover, FVC) , 537K M (Impervious sur- 
face area,ISA) 和 地 表 反 照 率 (Albedo)4 个 城市 地 表 
参数 用 于 研究 。 其 中 ,EVI 和 归 一 化 植被 指数 (Nor- 
malized difference vegetation index, NDVI) 是 当前 监 
测 地 表 植 被 变化 最 常用 的 2 种 植被 指数 ,大 量 研 究 
结果 表明 EVI 对 于 稀 叭 植被 的 监测 效果 优 于 NDVI， 
因而 更 适用 于 城市 区 域 植被 监测 “5 。 其 计算 公式 
如 下 : 


EVI=2.5 x 


(NIR - RED) 
(NIR + C, x RED - C, x BLUE Z) 

(10) 
式 中 :BLUE RED 和 NIR 为 Landsat 影 像 的 蓝光 、 红 
光 和 近 红 外 反射 率 波 段 ;L 为 土壤 背景 调节 参数 ;C， 
和 C0; 为 气 溶胶 阻抗 系数 。 

FVC 是 衡量 区 域 城市 地 表 植 被 状况 的 男 一 个 
重要 指标 " ,本 研究 拟 将 其 与 EVI 相 结合 一 同 用 于 
城市 夏季 地 表 温 度 变 化 分 析 。FVC 的 计算 公式 
如 下 : 


_[ (NDVI-NDVI) | 4 
awe ow T NDVI nin abd T ( 1 1 ) 
_ (NIR - RED) 
NDVI= NIR * RED) 


ISA (Ey Mili x E 38 AY RAMEN, de 
城市 地 表 热 环境 研究 的 关键 参数 。 研 究 选 用 3 Hg 
盖 研 究 区 的 夏季 末期 或 秋季 初期 Landsat 遥感 影像 ， 


获取 日 期 分 别 为 2001 年 9 月 12 日 .2011 年 10 月 3 日 
和 2018 年 9 月 11 日 ,以 此 避免 夏季 高 植被 覆盖 度 影 
响 城区 内 部 ISA 的 提取 精度 。 本 人 研究 基于 归 一 化 光 
谱 混 合 分 析 (Normalized spectral mixture analysis , 
NSMA) ”和 线性 光谱 混合 分 析 (Liner spectral mix- 
ture analysis,LSMA) “提取 了 研究 区 ISA 信 息 。 首 
先 , 利 用 NSMA 减少 不 同 用 地 类 型 端 元 的 光谱 差异 ， 
其 计算 公式 如 下 : 


=. AR N 
R,- 7100, 其 中 = 了 >R， (12) 
b=1 


式 中 : R, 表示 在 波段 8» 上 像 元 的 归 一 化 反射 率 ; u 
为 像 元 的 平均 反射 率 ; R, 表示 像 元 ! 的 原始 反射 
率 ;V 表 示 波 段 数 量 。 其 次 ,研究 利用 线性 光谱 解 混 
办 法 (LSMA) 计 算 城市 内 部 像 元 组 分 ,其 计算 公式 
如 下 : 


Se SN 
R= fxR,*e,, 其 中 > f=1;f<1 (13) 
i=] i=l 


式 中 : R,, KIED Be b 中 端 元 ;的 的 归 一 化 反射 率 
值 ; e, 表示 拟 合 误差 ; f 表示 端 元 i 的 权重 ;NN 表示 
端 元 数量 。 
反照 率 是 影响 城市 净 辐 射 以 及 能 量 平衡 的 一 
种 重要 因素 ,Mackey 等 研究 结果 表明 城市 地 表 
下 热 面 反射 率 越 低 其 LST 值 越 高 。 本 研究 运用 Li- 
ang ”基于 Landsat 卫 星 数据 提出 的 算法 提取 研究 区 
反照 率 ,其 计算 公式 如 下 : 
a =0.356x BLUE+0.130x RED+0.373 x NIR+ 
0.085 x SWIRI + 0.072 x SWIR2 - 0.018 
(14) 
式 中 : xu 表示 Landsat 数据 短波 反射 率 ; BLUE 
RED NIR,SWIRI1,SWIR2 分 别 表 示 Landsat 数 据 蓝 
光 、 红 光 、 近 红外 、 短 波 红外 1 波段 和 短波 红外 2 
波段 。 


3 结果 与 分 析 


3.1 城区 和 郊区 范围 的 界定 

SUHI 指 的 是 城区 与 郊区 温度 的 差 值 , 因 而 计 
算 研 究 区 夏季 SUHI 变 化 的 前 提 是 定义 研究 区 的 城 
区 范围 和 郊区 范围 。 本 研究 利用 2001、2011 年 和 
2018 年 获取 的 Landsat TM Fil Landsat 8 遥感 影像 ( 云 
履 善 度 <1% ) ,基于 支持 向 量 机 (Support vector ma- 
chine,SVM) 分 类 方法 ”提取 2001 一 2018 年 研究 区 
城区 和 郊区 的 范围 变化 信息 ,影像 获取 时 间 均 为 9 
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月 ,因为 这 一 时 期 研究 区 地 表 反 射 率 变化 较 其 它 时 


的 2018 年 土地 利用 分 类 信息 精度 ( 表 3)。 研 究 结果 


间 段 更 小 ,因而 更 有 利于 界定 城区 界限 。 首 先 将 乌 
鲁 木 齐 市 的 土地 利用 类 型 划分 为 植被 .水 体 、 裸 土 
和 建设 用 地 等 4 类 ,并 基于 SVM 分 类 方法 提取 了 研 
究 区 2001、2011 年 和 2018 年 的 土地 利用 信息 。 其 
次 ,利用 2018 年 的 谷歌 高 分 影像 作为 验证 数据 , 随 
机 选取 1400 个 验证 点 评价 基于 SVM 分 类 方法 提取 


表明 基于 SVM 提取 的 研究 区 土地 利用 信息 精度 能 
够 满足 后 续 人 研究 需要 。 

在 分 类 工作 完成 的 基础 上 ,乌鲁木齐 市 历年 城 
区 界限 均 以 高 密度 且 集 中 的 建设 用 地 为 划分 依据 ， 
郊区 界限 则 在 城区 界限 确定 后 分 别 围绕 城区 设 定 
一 个 较 小 的 缓冲 区 范围 (面积 占 市 区 面积 的 100%) 


表 3 基于 支持 向 量 机 (SVMD) 提 取 的 2018 年 土地 利用 信息 精度 


Tab.3 Accuracy of support vector machine classification in 2018 


验证 数据 


植被 水 体 HE ”建设 用 地 
植被 172 0 13 15 
水 体 0 20 0 0 
裸 土 7 1 135 37 
建设 用 地 30 0 37 933 
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图 2 乌鲁木齐 市 2001.2011 年 和 2018 年 城区 和 郊区 的 范围 


验证 数量 分 类 数量 正确 数量 ” 生产 者 精度 /% ”用户 精度 /% 
209 200 172 82.30 86.00 
21 20 20 95.24 100.00 
185 180 135 72.97 75.00 
985 1000 933 94.72 93.30 
87°20'E 87°40'E 


44°00'N 
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Fig.2 Urban and rural areas of Urumqi in the years 2001, 2011, and 2018 
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以 此 作为 2 个 不 同 郊区 范围 的 界限 ,从 而 开展 不 同 
郊区 界限 下 的 地 表 热 岛 强度 变化 监测 研究 。 如 图 2 
所 示 ,在 城区 和 郊区 范围 确定 后 ,最 终 确 定 了 一 个 
面积 为 50 kmx50 km 的 研究 区 范围 。 
3.2 Landsat LST 和 MODIS LST 数 据 的 融合 

人 研究 首先 通过 测试 实验 验证 FSDAF 模 型 用 于 
人 研究 区 LST 时 空 融合 的 实际 效果 。 图 3 所 展示 的 为 
用 于 时 空 融合 及 精度 验证 的 Landsat LST 和 
MODI11A1 ,以 及 基于 FSDAF 模型 融合 生成 的 “LST 
预测 数据 ”。 图 3a~b \i~j .na~o 为 用 于 模型 输入 的 基 
础 时 刻 Landsat LST fl MOD11A 1 影像 对 ,获取 时 间 
分 别 为 2001 年 5 月 23 日 .2011 年 6 月 4 日 和 2018 年 
5 月 22 日 ,图 3cfkpsv 为 预测 时 刻 的 MOD11A1 


输入 数据 ,以 此 融合 生成 2001 年 7 月 10 日 .2001 年 
9 月 12 日 .2011 年 8 月 7 日 .2018 年 6 月 23 日 .2018 
年 8 月 26 日 和 2018 年 9 月 11 日 高 空间 分 辨 率 的 
“LST 预 测 数据 "(图 3e hmsrx)。 研 究 采 用 预测 时 
刻 高 空间 分 辨 率 "LST 观 测 数据 “图 3dg lq tw) 
用 于 验证 “LST 预 测 数据 ”的 精度 。 

如 图 3 所 示 , 基 于 FSDAF 模 型 生成 的 “LST 预 测 
数据 ”与 同期 真实 获取 的 高 空间 分 辩 率 LST 观测 数 
据 ? 在 LST 总 体 空间 分 布 上 较为 一 致 。 从 表 4 可 以 
看 出 ,二 者 之 间 的 尼 ,AAD 和 RMSE 值 域 范围 分 别 在 
0.86~0.94 K ,2.15~2.59 K Fil 2.38~2.76 KK 之 间 。 

综 上 ,基于 FSDAF 模型 融合 生成 高 时 空 分 辩 率 
“LST 预 测 数据 ”的 精度 能 够 满足 后 续 研 究 需 要 。 
此 ,本 研究 同样 利用 测试 实验 中 用 于 模型 输入 的 基 


m LST/C m LST^C m LST/°C 
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(s) MOD11 


(t) Landsat 8 LST (u) LST 预测 数据 


(v) MODI1 (w)Landsat8 LST (x) LST 预测 数据 


3 “LST 观 测 数据 ”与 基于 FSDAF 模型 生成 的 “LST 预 测 数据 ” 


Fig. 3 Observed and predicted LST 


表 4 测试 实验 中 “LST 预测 数据 与 高 空间 分 辨 率 "LST 观 测 数据 "之 间 的 尼 ,AAD #1 RMSE f& 
Tab. 4 R’, AAD, RMSE values between the predicted LST and the actually observed Landsat LST of the testing experiment 


日 期 (年 -月 -日 ) 2001-6-10 2001-9-12 2011-8-7 2018-6-23 2018-8-26 2018-9-11 
R 0.93 0.88 0.94 0.91 0.90 0.86 
AAD/K 2.16 2.44 2.15 2.27 2.30 2.59 
RMSE/K 2.38 2.66 2.38 2.43 2.47 2.76 


注 :R, 决 定 系数 ;AAD ,绝对 平均 差 ;RMSE, 均 方 根 误差 。 通 过 0.01 水 平 显著 性 检验 (P<0.001)。 
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础 时 刻 Landsat LST 和 MOD11A1 影像 对 作为 FSDAF 
模型 基础 时 刻 (4) 的 输入 影像 ,将 8 d 合 成 的 夏季 
MOD11A2 数据 (获取 年 中 天 数 分 别 为 每 年 的 第 
169-176 d, 185-192 d, 193-1200 d, 201-1208 d, 
209-1216 d,217~1224 d,225-1232 d, 233-1240 d, 
241-1248 d All 249-1256 d) F HA TRU EN Za] Gn) 
的 MOD11A1 数 据 , 从 而 预测 生成 时 间 分 辨 率 为 8d 
覆盖 研究 区 的 高 空间 分 辩 率 夏季 LST 数 据 集 。 最 
后 ,通过 各 年 份 夏季 LST 数 据 之 和 求 平 均 ,获取 
2001.2011 年 和 2018 年 的 最 终 夏 季 LST 空 间 分 布 特 
征 ( 图 4)。 

如 图 4 所 示 , 乌 鲁 木 齐 市 区 面积 由 2001 年 的 
213.68 km 增加 至 2011 年 的 304.88 km’, 最 终 在 
2018 年 达到 了 533.33 km:。 随 着 城市 化 进程 加 快 ， 
特别 是 “ 乌 昌 一 体 化 ?战略 的 执行 ,城市 的 经 济 快速 
发 展 ,人 口 和 大 型 工业 企业 逐步 开始 向 城市 北部 迁 
移 。 因 而 在 2001 年 ,乌鲁木齐 市 夏季 LST 高 值 集中 
在 城市 中 心 区 域 。 而 在 2011 一 2018 年 ,夏季 LST 高 
值 则 主要 集中 在 城市 北部 ,并 且 逐 渐 向 东北 和 西北 
2 个 方向 扩张 。 

33 夏季 热岛 强度 和 LST 变 化 分 析 

研究 采用 SUHIL 和 SUHIL 等 2 个 热岛 强度 指标 
计算 研究 区 热岛 强度 (图 $)。 如 图 Sa 所 示 , 2001— 
2018 年 研究 区 的 热岛 强度 呈现 增加 的 趋势 。 其 中 ， 
SUHIL 由 2001 年 的 1.24 "C 增 加 到 2011 年 的 2.36 °C, 
2018 年 为 2.83 *C, SUHIL M) ri 2001 4Efff] 1.44 °C 1% 
Jil 8] 2011 4E BS 2.58 °C, 2018 4E 29 2.88 °C (图 5b)， 
表明 研究 区 郊区 范围 设 定 的 越 大 ,热岛 强度 就 越 
高 。 这 主要 是 由 于 乌鲁木齐 市 远郊 的 土地 利用 类 
型 以 植被 为 主 , 植 被 面积 占 比 随 着 郊区 范围 的 扩大 
而 增长 。 在 干旱 区 绿洲 城市 ,城区 内 部 植被 面积 的 
增加 有 助 于 缓解 城市 热岛 效应 。 而 在 城区 植被 条 
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图 4 2001 一 2018 年 研究 区 夏季 LST 时 空 分 布 特征 
Fig.4 Spatial distribution of predicted summer daytime 
LST during 2001—2018 


(b) 夏季 热岛 强度 均值 变化 


= SUHIL, 
= SUHII, 


2001 2011 2018 
年 份 


图 5 研究 区 2001.2011 年 和 2018 年 夏季 热岛 强度 (SUHITD) 变 化 
Fig. 5 Temporal changes of summer SUHII in the study area for the year 2001, 2011, and 2018 
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6 研究 区 2001 .2011 年 和 2018 年 夏季 不 同 土地 利用 lr SIUE 
类 型 的 地 表 温 度 (LST) 变 化 z z 
Fig. 6 Mean summer daytime LST by different land cover S P 
types in 2001, 2011, and 2018 + 
件 不 变 的 前 提 下 , 仅 郊区 植被 增加 则 会 导致 城市 地 
表 热 岛 强 度 增 加 。 
如 图 6 所 示 , 乌 鲁 木 齐 市 夏季 LST 值 在 不 同 土 z z 
© i) 
地 利用 类 型 上 的 变化 由 高 到 低 排序 依次 为 : 裸 土 > 中 
建设 用 地 > 植被 > 水 体 。 其 中 , 裸 土 LST 值 最 高 是 由 TE ee 
于 裸 地 的 土壤 水 分 含量 低 因 而 造成 热 惯量 低 。 相 87°20'E 87°40'E 
比 建设 用 地 , 裸 地 显 热 更 容易 消散 且 有 效 存储 的 能 _ 
量 较 少 ,因而 夏季 LST 值 最 高 。 植 被 的 LST 值 较 低 £ 3 
是 因为 植被 能 够 通过 冠 层 蒸腾 作用 有 效 释 热量 。 3 + 
水 体 由 于 具有 较 强 的 蓄 热 能 力 因 而 LST 值 最 低 。 在 "ms 
2001 年 和 2011 年 ,研究 区 建设 用 地 的 LST 值 随 着 郊 Bu 
区 范围 的 扩大 略 有 下 降 ,而 在 2018 年 则 呈现 增长 。 1 
这 表明 随 着 城市 化 进程 的 加 快 ,城市 远郊 的 乡镇 经 " z 
y H ~ x —- ` K © 
济 在 近 10 a 快速 发 展 ,远郊 乡镇 建设 用 地 面积 的 增 5 5 M, 
vt 
加 进一步 导致 了 研究 区 建设 用 地 LST 的 升 高 。 在 "o 87°208 87°40E 
(c) 20184 


2001 一 2018 年 ,研究 区 植被 的 LST 值 随 着 郊区 范围 
的 扩大 而 降低 , 裸 土 的 LST 值 同期 升 高 ,水 体 则 稳定 
不 变 , 这 表明 乌鲁木齐 市 远郊 的 土地 利用 类 型 以 植 
被 和 裸 土 为 主 。 因 此 ,在 干旱 区 绿洲 城市 城区 面积 
及 内 部 各 用 地 组 分 稳定 的 情况 下 , 当 城 市 郊区 的 土 
地 利用 类 型 以 植被 为 主 则 会 出 现 热岛 效应 ,而 当 郊 
区 土地 利用 类 型 以 裸 地 为 主 则 会 出 现 冷 岛 效 应 。 
3.4 夏季 LST 与 城市 地 表 参 数 之 间 的 定量 关系 
城市 地 表 参 数 变化 是 导致 干旱 区 绿洲 城市 夏 
季 LST 变 化 的 主要 因素 之 一 。 本 研究 基于 2001、 
2011 年 和 2018 年 夏季 Landsat 无 云 影像 分 别提 取 了 


图 7 研究 区 2001.2011 年 和 2018 年 夏季 增强 植被 指数 
(EVDT) 空 间 分 布 
Fig.7 Spatial distribution of EVI in 2001, 2011, and 2018 


EVI( 17) .FVC( I 8) F Albedo ( E] 9) ,基于 秋季 
Landsat 无 云 影像 提取 了 ISA( 图 10)。 

研究 首先 在 研究 区 城区 和 较 小 郊区 范围 内 随 
机 选取 1000 个 采样 点 。 其 次 ,在 较 大 郊区 范围 内 再 
随机 增 选 250 个 采样 点 (共计 1250 个 点 )。 最 后 , 利 
用 Pearson 相关 系数 分 析 研 究 乌鲁木齐 市 2001 一 
2018 年 夏季 LST 和 城市 地 表 参 数 的 定量 关系 ( 表 
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图 8 研究 区 2001.2011 年 和 2018 年 夏季 植被 覆盖 度 
(FVC) 空 间 分 布 
Fig. 8 Spatial distribution of FVC in 2001, 2011, and 2018 


5)。 如 表 5 所 示 , 乌鲁木齐 市 夏季 LST 和 EVI、FVC 
星 负 相关 ,表明 植被 的 增加 可 以 有 效 降 低 城市 地 表 
温度 ,改善 研究 区 整体 热 环 境 。 由 于 乌鲁木齐 市 的 
东部 为 林地 ,南部 为 草地 ,北部 多 为 耕地 ,因而 随 着 
城区 的 扩张 和 郊区 范围 的 扩大 LST 与 植被 之 间 负 相 
关 性 更 为 显著 。 另 一 方面 ,乌鲁木齐 市 夏季 LST 和 
ISA 呈正 相关 ,在 2001 年 和 2011 年 二 者 之 间 的 Pear- 
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图 9 研究 区 2001,.2011 年 和 2018 年 夏季 地 表 反照 率 
(Albedo) 空 间 分 布 
Fig.9 Spatial distribution ofAlbedo in 2001, 2011, and 2018 


son 相关 系数 随 着 郊区 范围 的 扩大 而 降低 ,在 2018 
年 则 呈现 增加 的 趋势 ,这 主要 是 由 于 近年 来 城市 远 
郊 乡 镇 城市 化 进程 加 快 导致 不 透水 面 面 积 呈 快速 
增加 。 最 后 ,由 于 人 研究 区 不 同 郊区 范围 内 均 有 大 量 
裸 土 , 干 旱 区 裸 土 仿 水量 低 且 反照 率 高 ,因而 与 
LST 呈 正 相 关 , 且 相关 性 随 着 郊区 范围 的 扩大 而 
增强 。 
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图 10 研究 区 2001、2011 年 和 2018 年 夏季 不 透水 而 (ISA ) 空 间 分 布 
Fig. 10 Spatial distribution of ISA in 2001, 2011, and 2018 


表 5 2001,.2011 年 和 2018 年 研究 区 夏季 LST 和 EVIFVC ISA, Albedo 的 Pearson 相关 系数 
Tab.5 Pearson correlation coefficients between summer daytime LST and EVI, FVC, IS, and Albedo in 2001, 2011 and 2018 


2001 4F 201146 2018 年 
Area 1 Area 2 Area l Area 2 Area 1 Area 2 
EVI -0.79 -0.81 -0.81 -0.82 -0.84 -0.85 
FVC -0.67 —0.68 -0.73 -0.77 -0.74 -0.80 
ISA 0.85 0.82 0.88 0.86 0.88 0.89 
albedo 0.86 0.88 0.86 0.89 0.91 0.92 


注 :EVI 为 增强 植被 指数 ;FVC 为 植被 覆盖 度 ;ISA 为 不 透水 面 ;Albedo 为 地 表 反 照 率 ;Area 1 表示 城区 和 较 小 郊区 范围 ;Area 2 表示 城区 和 较 
大 郊区 范围 。 通 过 0.01 水 平 显 著 性 检验 (P<0.001)。 


增加 的 趋势 , 且 郊 区 范围 越 大 ,热岛 效应 越 明 显 ; 

4 结论 (3) 在 各 土地 利用 类 型 中 , 裸 地 的 夏季 IST 最 高 ,其 
本 研究 以 于 时 区 典型 绿洲 城市 乌鲁木齐 次 为 建设 用 地 和 植被 ,水体 最 低 ; (4) 研究 区 夏季 
市 为 研究 区 ,基于 FSDAF 模 型 生成 了 2001 .2011 年 。 ST 升 高 与 地 表 反 照 率 和 1SA 的 增加 呈正 相关 ,与 
和 2018 年 夏季 高 时 空 分 状 率 LST 数 据 集 ,并 据 此 开 ”植被 因子 呈 负 相关 ;(5) 在 二 时 区 绿洲 城市 ,城区 内 
展 了 乌鲁木齐 市 夏 委 地 表 热 岛 效 应 变化 研究 。 研 。 部 植被 面积 的 增加 有 助 于 缓解 城市 热岛 效应 ,而 仪 
究 结果 表明 :(1) 基于 FSDAF 模 型 融合 生成 的 "LST ” 郊区 植被 的 增加 则 会 导致 城市 地 表 热 岛 强 度 增 
预测 数据 ”精度 较 高 ,能够 满足 研究 需要 ;(2) 乌 鲁 。 加 。 因 此 ,当地 政府 绿化 部 门 应 在 保证 郊区 植被 数 
木 齐 市 夏季 热岛 强度 在 不 同 的 郊区 范围 内 均 呈 现 量 的 前 提 下 ,不 断 提 升 城市 内 部 的 植被 种 植 数量 ， 
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提高 植被 覆盖 度 ,以 此 减缓 城市 热岛 效应 的 发 生 
发 展 。 
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FHEWY 


Spatiotemporal variations of the summer daytime surface urban heat island 
of oasis city in arid area based on FSDAF model 


WANG Shuang, WANG Chengwu, ZHANG Feiyun 
(School of Management, Xinjiang Agricultural University, Urumqi 830052, Xinjiang, China) 


Abstract: The urban heat island effect and ecological and environmental problems of an oasis city in an arid area 
resulting from the acceleration of the global urbanization process have become one of the current research 
hotspots worldwide. As a critical parameter in investigating the surface urban heat island (SUHI), the land surface 
temperature (LST) retrieved from thermal infrared images is, to some extent, inaccurate because there is no single 
sensor that can capture real-time LST at both high spatial and temporal resolutions. This paper used the flexible 
spatiotemporal data fusion method to generate high spatiotemporal resolution summer daytime LST data in 2001, 
2011, and 2018. Spatiotemporal variations of the summer daytime SUHI over Urumqi, an oasis city in the arid 
area of western China, were assessed based on several urban surface biophysical variables. The results show that 
the SUHI intensity (SUHII) in Urumqi during the study period was calculated using two indicators of SUHII, (the 
urban and smaller rural area difference in the average LST) and SUHIL (the urban and larger rural area difference 
in the average LST). Significantly increasing trends of SUHII in the study area were observed. SUHII, increased 
from 1.24 °C in 2001 to 2.36 °C in 2011 and 2.83 °C in 2018, whereas SUHIL increased from 1.44 °C to 2.58 °C 
and 2.88 °C in the same periods. The highest and lowest summer daytime LST values were observed over areas 
of bare soil and water, respectively. The distribution of the summer daytime LST correlated positively with the 
albedo, the impervious surface area, and it correlated negatively with the enhanced vegetation index and 
fractional vegetation cover. The results emphasize the role of bare soils in aggravating the SUHI in cities in arid 
areas. Finally, we find that in oasis cities in arid areas, such as Urumqi, although increasing the amount of 
vegetation covering the urban area may be an effective way to mitigate the SUHI, more profuse vegetation 
coverage within a larger rural area will increase the SUHII during the summer daytime. 


Key words: surface urban heat island; spatiotemporal fusion model; oasis city in arid area 


